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'-^ABSTRACT 

1-H Trajectory planning is a critical step while programming the 
^ parallel manipulators in a robotic cell. The main problem arises 
when there exists a singular configuration between the two poses 

QQ of the end-effectors while discretizing the path with a classical 
approach. This paper presents an algebraic method to check 
the feasibility of any given trajectories in the workspace. The 
solutions of the polynomial equations associated with the tra- 
jectories are projected in the joint space using Grobner based 
elimination methods and the remaining equations are expressed 
in a parametric form where the articular variables are functions 
^ of time t unlike any numerical or discretization method. These 
formal computations allow to write the Jacobian of the manip- 
ulator as a function of time and to check if its determinant can 
vanish between two poses. Another benefit of this approach is 
to use a largest workspace with a more complex shape than a 
cube, cylinder or sphere. For the Orthoglide, a three degrees 
of freedom parallel robot, three different trajectories are used to 
illustrate this method. 


INTRODUCTION 

One of the crucial steps in the trajectory planning is to check 
the singularity-free paths in the workspace for the parallel ma¬ 
nipulators. It becomes a necessary protocol to validate the tra¬ 
jectory when the parallel robot is used in practical applications 


such as precise manufacturing operations. A trajectory verifi¬ 
cation problem is presented in Q based on some validity crite¬ 
ria like whether the trajectory lies fully inside the workspace of 
the robot and is singularity-free. Singularity-free path planning 
for the Gough-Stewart platform with a given starting pose and a 
given ending pose has been addressed in using the clustering 
algorithm is presented in . An exact method and an approxi¬ 
mate method are described in 0 to restructure a path close to the 
singularity locus into a path that avoids it while remaining close 
to the original path. Due to the geometry of the mechanism, the 
workspace may not cover fully the space of poses 0, hence it is 
necessary to analyze the workspace of the manipulator. 


The workspace of a parallel robot mainly depends on the ac¬ 
tuated joint variables, the range of motion of the joints and the 
mechanical interferences between the bodies of the mechanism. 
There are different techniques based on geometric tools 00, 
discretization ||7}0, and algebraic methods |T0}fT^ which are 
used to compute the workspace of a parallel robot. An algebraic 
method to solve the forward kinematics problem specifically ap¬ 
plied to spatial parallel manipulators is described in |T4| . The 
main advantage of the geometric approach is that, it establishes 
the nature of the boundary of the workspace (T5). A procedure 
to automatically generate the kinematic model of parallel mech¬ 
anisms which further used for singularity free path planning is 
reported in GD- 

An algorithm for computing singularity-free paths on 
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closed-chain manipulators is presented in p7| , also this method 
attempts to connect the given two non-singular configurations 
through a path that maintains a minimum clearance with respect 
to the singularity locus at all points. The main drawback of any 
numerical or discretization methods is that there might be a sin¬ 
gular configuration between two poses of the end-effectors while 
discretizing the path. This paper illustrates a technique based on 
some algebraic methods to check the feasibility of any given tra¬ 
jectories in the workspace : it allows to write the Jacobian of the 
manipulator as a function of time and to check whether its deter¬ 
minant vanishes between two poses. Also, when the trajectory 
meets a singularity, its location can also be computed. 

The outline of this paper is as follows. We first introduce the 
modelization and the basic algebraic tools for computing trajec¬ 
tories in the workspace for parallel robots. We then extended the 
process to check the singularity-free paths within the workspace 
and ensures the existence of a singular configurations between 
the two poses of the end-effector by analysing the Jacobian as a 
function of time. In later sections we give some detailed exam¬ 
ples using three different trajectories for the Orthoglide, a three 
degrees of freedom parallel robot. 


METHODOLOGY 

To ensure the non existence of a singular configurations be¬ 
tween two poses of the end- effector, it is necessary to express 
the Jacobian of the manipulator as a function of the time. The 
general procedure to check the feasibility of the trajectories in 
the workspace as well as a function to compute the projection of 
the trajectories in the joint space can be decomposed as follows : 

• Defining the constraint equations, articular and pose vari¬ 
ables associated with the parallel mechanism; 

• Computing the singularities and their projections in the 
workspace and in the joint space; 

• Computing the workspace and joint space boundaries; 

• Computing a parametric form of the trajectories in the Carte¬ 
sian space; 

• Computing the projections of the trajectories in the joint 
space; 

• Computing the Jacobian of the manipulator as a function of 
the time. 

By computing, we mean, by default, getting a full charac¬ 
terization as solutions of an exact system of algebraic equations. 
By abuse of language, computing a projection might thus mean 
computing the algebraic closure of the projection, for example 
by using classical elimination strategies based on Grobner bases 
(which thus could introduce a subset of null measure made of 
spurious points). 

Kinematics involves the study of the position, velocity, ac¬ 
celeration, and all higher order derivatives of the position vari¬ 
ables (with respect to time or any other parameter/variables). 


Hence, the study of the kinematics of manipulators refers to all 
the geometrical and time-based properties of the motion. For a 
translational manipulator, the set of (algebraic) relations that con¬ 
nects the input values p to the output values X will be denoted 
by 


F(p,X)=0 


( 1 ) 


where p is the set of all actuated joint variables and X is the set 
of all pose variables of the end-effector. 

For the study of manipulators, two problems must be con¬ 
sidered: the direct kinematic problem (DKP) and the inverse 
kinematic problem (IKP). Specifically, given a set of joint val¬ 
ues, solving the DKP consists in computing the position and the 
orientation of the end-effector relatively to the base, or, equiv¬ 
alently to change the representation of the manipulator position 
from a joint space description into a Cartesian space description. 
Given the position and orientation of the end-effector of the ma¬ 
nipulator, solving the IKP problem consists in computing all the 
possible sets of joint angles or parameters that could be used to 
attain this position and orientation, or, equivalently, mapping lo¬ 
cations in three-dimensional Cartesian space to locations in the 
robot’s internal joint space. 

When considering an algebraic modelization, the direct 
kinematics has basically several solutions, which refers to the 
several poses of the end-effector for given values of the joint co¬ 
ordinates. It is therefore possible to assemble the manipulator 
in different ways, and these different configurations are known 
as the assembly modes of the manipulator |T^ . Similarly, the 
multiple inverse kinematic solutions induce multiple postures 
for each leg of the manipulator and is termed as the working 
modes of the manipulator. An analytical relation can be given as 
p = 7(X) for IKP whereas X = j8(p) for DKP. Differentiating 
Eq.0 with respect to time leads to the velocity model: 


AX + Bp =0 


where A = , B = — 

oX op 


( 2 ) 


The matrices A and B are respectively the direct-kinematics 
and the inverse-kinematics Jacobian matrices of the manipulator. 
These matrices are used for characterizing different kinds of sin¬ 
gularities. The parallel singularities occur whenever det(A) = 0 
and the serial singularities occur whenever det(B) = 0. The par¬ 
allel singular configurations are located inside the workspace. 
They are particularly undesirable because the manipulator can¬ 
not resist to any forces and its control is lost. 

Eliminating p in the system F(p,X) =0, det{A) = 0 by 
means of a Grobner basis computation for a suitable elimination 
ordering (see p9} ) defines (the algebraic closure of) the projec¬ 
tion ^ (X) of the parallel singularities in the workspace. In the 


same way, one can compute (the algebraic closure of) the projec¬ 
tion e(X) of the parallel singularities in the joint space. Both are 
then defined as the zero set of some system of algebraic equations 
and we assume that the considered robots are generic enough so 
that both are hypersurfaces. 

The set of equations associated with the joint limits of the 
actuator x(p) can also be projected, with the same elimination 
method, in the workspace /i(X) as in Eq. where pmin, pmax 
are the minimum and maximum values of the articular variables: 
X{p) and /i(X) are the crucial parameters in determining the 
number of assembly modes and working modes of the parallel 
manipulator. 


X{P) ^ ^P ^ (Pmin^Pmax] (3) 

To use algebraic tools such as Grobner bases, we must rep¬ 
resent the trajectories in an algebraic form. A classical approach 
can be used to transform the trigonometric equations to algebraic 
ones as in | [2Q| . Equation Q represents the trajectory in paramet¬ 
ric form within the workspace as the function of time t. 


X^0(t) (4) 

'T(X,p,t) in Eq. ^ is the system of equations which con¬ 
tains kinematic equations E(p,X) and the parametric equations 
of trajectory 0(t). 


'T(X,p,t) = [0(t)-X,F(p,X)] (5) 

This system of equations is projected in the joint space as 
T(p,t) as shown in Eq. ^ using Grobner based elimination 
method. 


'T(X,p,t)^T(p,t) (6) 

The parametric equations of the trajectory as a function of 
the time in the joint space can then be obtained by solving p ^ 
T(p,t). 

The workspace analysis allows the characterization of the 
workspace regions where the number of real solutions for the 
inverse kinematics is constant. A cylindrical algebraic decompo¬ 
sition (CAD) algorithm is used to compute the workspace of the 
robot in the projection space X with joint constraints x(p) taken 
in account | |2TH23l . 

The three main steps involved in the analysis are: 


• Computation of a subset of the joint space (resp. workspace) 
where the number of solutions changes: the Discriminant 
Variety . 

• Description of the complementary of the discriminant va¬ 
riety in connected cells: the Generic Cylindrical Algebraic 
Decomposition (CAD). 

• Connection of the cells belonging to the same connected 
component in the counterpart of the discriminant variety: in¬ 
terval comparisons. 

The joint space analysis predicts the feasible and non- 
feasible combinations of the prismatic joint variables which are 
essential for the parallel robot control. The joint space analysis 
allows the characterization of the regions where the number of 
real solutions for the direct kinematic problem is constant. The 
joint space analysis is done using CAD, which gives the num¬ 
ber of cells corresponding to the number of solutions in the joint 
space. x{p) parameter significantly changes the number of as¬ 
sembly modes and working modes of the manipulator. 

Substituting (j) (t) in ^ (X) (in the polynomial defining ^ (X) 
) we then get (0(t)), which defines the vanishing values of the 
Jacobian as a function of the time. This equation in cos(t),sin(t) 
can be turned into a zero-dimensional bivariate system and its 
solutions can thus be computed exactly, either by means of a 
rational parametrization or by means of isolating intervals with 
rational bounds that can be refined to any arbitrary precision 
Whatever the chosen (exact) representation for a solution, 
it can easily be checked if it vanishes between the two poses of 
the end-effector. The solutions of (0(t))contain the singular 
points on the trajectory 0(t) and eventually few spurious points 
due to the projection of det(A) in Cartesian space. Spurious sin¬ 
gular points can then be differentiated from real singular points 
by substituting 0(t) and the solutions of T(p,t) in det(A). 

METHOD VALIDATION 

We propose to validate our approach (and related tools) by 
checking the feasibility of three different trajectories for the 0th- 
oglide. The first two trajectories are heart shaped planar paramet¬ 
ric curves and the third trajectory is a parametric curve in three 
dimensions. 

These trajectories are fully inside the workspace and are se¬ 
lected such that their projections in the joint space is defined by 
parametric equations containing trigonometric functions. In ad¬ 
dition, we have shown that the singular and singular-free trajec¬ 
tories are well discriminated, as one of the trajectory cuts the 
singularity surface. 


^for example the functions Groebner[RationalUnivariateRepresentation] and 
RootFinding [Isolate] in Maple software 




FIGURE 1. Architecture of an Orthoglide including two assembly 
modes (a) Workspace plot with x(p) joint constraints (b) and Paral¬ 
lel singularity det(A) projected in the workspace as ^ (X) (c) and joint 
space plot with x{p) joint limits (d) 



FIGURE 2. A comparison between the parallel singularity surface for 
an Orthoglide computed with the joint limits jU(X) (a) without jU(X) (b) 


where I is the leg length of the manipulator and is equal to two 
for all the computations. The set of equations associated with 
the joint limits of the actuator x{p)^ projected in the workspace 
/i(X), are given in Eq. §. X{p) plays an important role in de¬ 
termining the shape of the workspace and singularity surfaces. It 
also affects the number of solutions for the IKP Le. the work¬ 
ing modes associated with the manipulator. Figure [^b) repre¬ 
sents the workspace of the Orthoglide. A CAD algorithm is used 
to compute the workspace of the robot in the projection space 
(x,y,z), taking into account the joint constraints z(p) • Without 
considering the joint limits, the Orthoglide admits two assembly 
modes and eight working modes p4| . 


Manipulator Architecture and Kinematics 

The manipulator under the study is an Orthoglide parallel 
robot with three degrees of freedom. The mechanism is driven by 
three actuated orthogonal prismatic joints pi, made of three par¬ 
allelogram connected by revolute joints to the tool center point 
on one side and to the corresponding prismatic joint at another 
side. The assembly modes of these robots depends on the solu¬ 
tions of the DKP as shown in Fig.[2a). The point P/ represents 
the pose of corresponding robot. However more than one value 
of i for the point denote multiple solutions for the DKP. 
is equal to pi, where pi represents the prismatic joint variables 
whereas P represents the position vector of the tool center point. 

The constraint equations for the Orthoglide with the actuated 
variables p = [pi,p 2 ,P 3 ] and the pose variables X = [v,y,z] are: 


F(p,X):(x-Pi)2+/+z2 = /2 
x^^{y-p2)^+z^ = l^ 
v2+/ + (z-P3)" = /" (7) 


x(p) = [pi,p2,p3,-4 + pi,-4 + p2,-4-f P3] 

X{p)^p{X) 

P; i-A — 4 /= 1,2,3 

-4 + pi - 8v-i-12 

—4 Pi i-A — 8y + 12 

—4 + p/i-A +y^ + z^ — 8z + 12 (8) 

In the Figure [^b), the yellow region (resp. green region), 
corresponds to the region where the inverse kinematic model has 
two real solutions (resp. height). 

The joint space analysis is done using CAD to hnd the re¬ 
gions where the direct kinematics problem admits a hxed number 
of real solutions. For example, in Fig. ^d), the red cell corre¬ 
sponds to to the region where the DKP has two solutions. 

There is a difference in the shape of the singularity sur¬ 
face, depending if we consider the joint constraints p(X) or not 
(Fig.|g. 























det(A) = - 8 piP 2 P 3 + 8 piP 2 Z + 8 piP 3 >> + 8 p 2 P 3 X 
det(A) !->• e(p) det(A) !->• (^ (X) 
e(p) = Pfp| + Pfpf + P^P^ + 2p?p|p| + P?P3' + 
P2V3 + P2P3 - 16pfp| - 16pfp| - 16p|p| 


(9) 


In Eq. det(A) is the parallel singularity of the Orthoglide 

and ^ (X) is the projection of det(A) in workspace. The mathe¬ 
matical expression for ^ (X) is not displayed in Eq. El due to 
the lack of space. Eig. [^c) shows the projection ^ (X) which is 
plotted with fi (X) as one of the input parameter. The degree of 
this characteristic surface is 18 and it represents the singularities 
associated with the eight working modes. 

Trajectory definition in the workspace 

Trajectories 1 & 2 are heart shaped parametric curves which 
( Eig. I^a)). Eq. ( p^ and Eq. ( pT] ) are the mathematical defi¬ 
nitions of Trajectory 1 and Trajectory 2, respectively. As these 
equations are trigonometric equations, it is necessary to represent 
them in an algebraic form. 0 i(r) and 02(0 the parametric 
equations for Trajectory 1 and Trajectory 2, respectively, which 
are defined for r G [—;r, ;r]. Erom now, sin and cos are replaced 
by s and c respectively to reduce the size of the expressions. 


01 (0 : 


8 3. s 

■= 


z = 1 


( 10 ) 



FIGURE 3. Position of Trajectories 1&2 (a) Trajectory 3 (b) in the 
workspace of an Orthoglide 


In order to turn the system to an algebraic one, we add the 
following equations 


sm{t) = sin J cos{t) = cosJ 
sin + cos = 1 


(13) 


and we remark that this change of variables does not introduce 
any spurious solutions since it is bjective (t e [—71, n]). 

Projection in the joint space 

01 (O’ ^ 2(0 ^ 3(0 trajectories which are defined 

in the workspace. To project these trajectories in the joint space, 
it is necessary to formulate a system of equations corresponding 
to each trajectory which also consists the kinematic equations of 
the manipulator. 'Pi, ^^2 and 'Ps are the corresponding systems 
of equations for the Trajectories 1, 2 and 3, respectively (see Eq. 

(GD¬ 


IS 


4 , 

(kit): x=-s^it) 


'Pi(x,p,t) = [<;)i(t)-x,F(p,x)] 


J 

1 1 1 


'P2(X,p,t) = [<fc(t)-X,F(p,X)] 


-c(2t) - c(3t) - c(4t) 

4 ^ ^ 20 ^ ^ 20 ^ ^ 


'F3(X,p,t) = [(fe(t)-X,F(p,X)] 

(14) 

z= 1 

( 11 ) 




Trajectory 3 is a parametric curve in three dimensions. 
Eig. I^b) represents the trajectory in the workspace of the ma¬ 
nipulator. The parametric equations of the trajectory is defined 
by Eq. (H^. is the parametric equation for Trajectory 3 and 
is defined for r G [ 0 , 20 ]. 


piit): x = s(t) y = c{t) z=—t ( 12 ) 


Each system of is projected in the joint space as Ti(p,t), 
T 2 (p,t) and T 3 (p,t) (see Eq. (p3])). By solving Ti(p), T 2 (p) 
and T 3 (p) for p, we get the corresponding parametric equations 
for the trajectories in the joint space, as shown in Eig. for the 
Trajectories 1&2 and in Eig. for the Trajectory 3. Eqs. G^- 
( p^ are used to plot all possible images of the trajectories in the 
joint space. All the computations are done without considering 
the joint limits. As the Orthoglide has eight working mode, there 
exists eight possible trajectories in the joint space ( Eig. and 
Eig.for Trajectories 1&2 and Trajectory 3, respectively). 












FIGURE 4. A pictorial representation of the mapping of trajectories from workspace to joint space. Eight different pairs of trajectories, Ti (p, t) & 
T 2 (p,t) in joint space are the image of corresponding 0i(t) & 02(t)- These eight different trajectories are associated with the eight working modes of 
the Orthoglide. Only one trajectory lies inside the joint space boundary due to to the joint constraints. 


'Pi(X,p,t)H^Ti(p,t) 

Vf e [-n, k] 


'P2(X,p,t)H^T2(p,t) 

Vf e [-TT, n] 


'I'3(X,p,t)^T3(p,t) 

Vf G [0, 20] 

(15) 

Solving Ti(p,t) andT 2 (p,t) 

for p gives eight possible so- 


lutions ( Eq. ( p^ for Trajectory 1 and Eq. ( p^ for Trajectory 2), 
which inferred eight different possible images of Trajectory 1&2 
in joint space. But from Eig.j^one can see that only one trajec¬ 
tory lies inside the joint space boundary for the both trajectories 
in workspace. 

Pi = s(r) ± 2 ^400 s2(/) - /2 +1200 

P 2 = c(f) ±2^400c2(r)-r2 + i200 

p3 = —t±V3 (16) 



FIGURE 5. Mapping of the Trajectory 3 from workspace to joint 
space. Eight different possible solutions of T 3 (p,t) are marked in joint 
space which are the image of 03 (t). There is only one feasible trajectory 
(marked as 1) lies inside the joint space boundary. 


Eigurej^ shows all the possible images of Trajectory 3 in the 
joint space. Solving T 3 (p,t) for p gives eight possible solutions 
which (Eq. (p^). It can be seen in the Eig.j^that only one trajec¬ 
tory (which is marked as 1), lies inside the joint space boundary. 


Singularity analysis 

To check if there exists any singular configuration between 
the two poses of the end-effector it is necessary to express the 
Jacobian of the manipulator as a function of the time or of some 






















Pi = ^s^(r) ± ^ ^1996s6(r) - 400s8(r) + 560s6(f)c(0 - 100s4(f)c(f) - 2009s4(f) + 2190c3(f) - 1379c2(r) - 1320c(f) + 3988 
P 2 = - ^c4(r) - ^c\t) - + gc(r) + ? ± ^ ^6Ac^{t)-\92c\t) + \92c\t) + %3 

P 3 = 1 ± ^^3200-400c8(r) -560 c 7(?) + 1204c6(/) + 1580c5(f) -3221c4(r) +710c3(/) + 3051c2(/) - 860c(0 (17) 


Pi = ^s^(r) ± ^^76s6(f) - 16s\t) + I6s^{t)c{t) + Us^{t)c{t)-S5s^{t)+ 96c^{t)-66c^{t)-60c{t)+ 321 
Pi = -lc\t) - - lc 2 (r) + ^c (0 + \ ±l^l 6 c 6 (r)- 48 c 4 ( 0 + 48 c 2(0 + 59 

P 3 = 1 ± :j^^332 - 16c8(r) - 16c7(f)’ + 52c6(r) + 60c5(f) - 145c4(r) + 24c3(f) + 132c2(r) - 32c(f) (18) 


independent variable. The proposed algebraic method enables us 
to write the Jacobian of the manipulator as a function of the time 
and to check if its determinant vanishes between two poses. 


By substituting the values of 0i(^), 02(0 03(0 from 

Eq. Eq. ([H]) and Eq. ( p^ in /i(X), we will get /i(0i(O), 
M (02 (0) M (03 (0) the determinant of the Jacobian for Tra¬ 

jectory 1, Trajectory 2 and Trajectory 3 respectively (Eq. <lig). 
Due to the large expressions, the equations for /ii(t), /i2(t) and 
/i 3 (t) are not presented but their roots define the singular configu¬ 
rations : Eiguresj^ [T^and pT]show the values of these functions 
when t varies. 


pi(t) = p((/»i(t)) 

Piit) = Pi<hit)) Vte[-;T,7r] 

P3(t) = M(<^)3(t)) Vte[0, 20] 


(19) 



FIGURE 6. Representation of the trajectories with singularity sur¬ 
face. Trajectory 1 cuts the parallel singularity surface <^(X) in four 
points si, S2, S2 and ^4 in the workspace of the Orthoglide. Also, it can 
be seen that Trajectory 2 lies inside ^ (X) as 1^2 W Vt G [—tt, tt]. 


The number of solutions of Eq. ( p^ gives the total num¬ 
ber of singular points on the corresponding trajectory. The so¬ 
lutions for /ii(t), /i2(t) and /i 3 (t) are shown in Eq. ( [^ . Erom 
Eq. (20) we get that there exists four solutions for /ii(t) and 
zero solutions for /i2(t) and /i 3 (t), which confirms the presence 
of singular points on Trajectory 1 whereas Trajectory 2&3 are 
singularity-free trajectories. Note that numerical approximations 
are given for readability, but, in practice, isolating intervals with 
rational bounds with arbitrary width can be computed (see sec¬ 
tion METHODOLOGY or |25 chapter 8] for more details) so 


that the result is certified. 


t= [-1.51,-0.97,0.97,1.51] ^ All (t) =0 
t= [0.97,1.51] ^det(A) =0 

t = [ ] ^ /i2 (t) = 0 

t = [ ] ^ M3 (t) = 0 


Vt G [-71,71] 

Vt G [-71,71] 

Vt G [-71,71] 

Vt G [0, 20](20) 


/ii(t) is the determinant of the Jacobian corresponding to 
the Trajectory 1. By solving a zero-dimensional system of equa- 




















FIGURE 7. Representation of the real singular points 53 and >^4 along 
with the singularity surface ^ (X), which is associated with one working 
mode. 


tions, it can be shown that Trajectory 1 cuts the singularity sur¬ 
face in four points S 2 , and S 4 (see Fig.j^and Eq. ([20|)). The 
image of these points in workspace is computed by substituting 
the values of Eq. ( [20| ) in Eq. Similarly, the image of these 
points in joint space is computed by substituting the values of Eq. 
( [^ in Eq. ([T7](. The obtained values are tabulated in Table [T] 
Eigures represents the Trajectory 1&2 along with singularity 
surface <^(X), also, all the singular points Si, S 2 , S 3 and S 4 are 
located on the singularity surface. 

Si and S 2 are the spurious singular points which were intro¬ 
duced by the projection of det(A) in Cartesian space. Substitut¬ 
ing /ii(t) and 01 (t) in det(A) gives two solutions ( Eq. (20)). 
These two solutions S 3 and S 4 are the real singular points on Tra¬ 
jectory 1. The variation of det(A) along Trajectory 1 is shown in 
Fig.0 



FIGURE 8. Variation of jUi(t) along Trajectory 1. ^i, ^2, ^3 and >^4 
represents the four solutions of fii (t) = 0 Vt G [—tt, tt]. 



FIGURE 9. Variation of det(A) along the Trajectory 1. 53 and 54 
represents the two solutions of det(A) = 0 Vt G [—tt, tt]. 


TABLE 1 . Singular postures on Trajectory 1 


Workspace 

Joint space 

S 

X 

3^ 

z 

Pi 

P 2 

P 3 

Si 

-1.13 

0.35 

1.00 

0.55 

1.66 

2.60 

Si 

-0.65 

0.80 

1.00 

0.88 

2.40 

2.71 

S3 

0.65 

0.80 

1.00 

2.18 

2.40 

2.71 

S4 

1.13 

0.35 

1.00 

2.83 

1.66 

2.60 


The true singular postures S 3 and S 4 , belonging to the singu¬ 
larity surface <^(X), which is associated with the working mode 
used by the Orthoglide prototype ( Eig.j^. /i2(t) and /i 3 (t) are 
the determinanst of the Jacobians corresponding to the Trajecto¬ 
ries 2&3 respectively. 

The variation of /i2(t) along Trajectory 2 for r G [— TT, ;r] is 
shown in Fig.[^ Figure pT] represents the variation of /i 3 (t) for 
t G [0,20] along Trajectory 3. All the values given in Tableare 
associated with the trajectories which is marked as 1 in Fig. 
These values are obtained by solving a zero-dimensional system 
of polynomial equations using RootFinding:-Isolate so that the 



























FIGURE 10 . Variation of jU2(t) along the Trajectory 2 . There are no 
solutions for jU2(t) =0 Vt G [— tt, tt]. 


Conclusions 

This paper reports the use of algebraic methods to check the 
feasibility of given trajectories in the workspace. The general 
method uses the projection of the polynomial equations associ¬ 
ated with the trajectories in the joint space using a Grobner based 
elimination strategy. There is a significant change in the shape 
of workspace and the singularity surfaces due to the joint con¬ 
straints. The proposed method enables us to project the joint con¬ 
straints in the workspace of the manipulator, which further helps 
to analyse the projection of the singularities in the workspace 
with the joint constraints. Such computations ensure the exis¬ 
tence of a singular configuration between two poses of the end- 
effector unlike other numerical or discretization methods. This 
paper highlights that the singularity analysis should be done in 
the cross product of the workspace and joint space for parallel 
robots with several assembly and working modes. The single 
analysis of the singularities in a projection space can introduce 
spurious singularities. In fact, the algebraic tools does not allow 
to distinguish the parallel singularities according to the working 
mode. 



FIGURE 11 . Variation of jU3(t) along the Trajectory 3 . There are no 
solutions for =0 Vt G [—tt, tt]. 


related values for sinJ and cosJ are certified. 

In Table Si and S 2 are the spurious singular points and S 3 
and S 4 are the real singular points on the Trajectory 1. Trajectory 
2&3 are the singularity-free trajectories. 
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